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Abstract 

Results on hadronic final states in deeply inelastic scattering are reviewed. They comprise jet 
production and its interpretation in perturbative QCD, signatures to distinguish conventional 
QCD dynamics from possible new features of QCD at small x, and measurements of inclusive 
charged particle production. Theoretical developments such as color dipole emission and instanton 
induced final states are reported on. 



1. Introduction 

The basic measurement in deeply inelastic scattering 
(DIS) is a measurement of the cross section ep — > eH in 
terms of the structure function F 2 , where H stands for 
any hadronic system. A wealth of information upon the 
partonic structure of the proton and its dynamics have 
been obtained from structure function measurements. 
Measurements of the properties of the hadronic final 
state H provide complementary information which 
cannot be obtained from inclusive structure functions. 

In the simple quark parton model (QPM) of DIS, a 
quark is scattered out of the proton by the virtual boson 

* Summary talk from the working group II "Hadronic Final 
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emitted from the scattering lepton. QCD modifies this 
picture. Partons may be radiated before and after the 
boson-quark vertex, and the boson may also fuse with a 
gluon inside the proton by producing a quark-antiquark 
pair (figure |l|). In fact, the parton which is probed by 
the boson may be the end point in a whole cascade of 
parton branchings. This parton shower materializes in 
the hadronic final state, allowing experimental access to 
the dynamics governing the cascade. 

HERA has opened a new kinematic domain to study 
QCD in DIS, and most contributions in this working 
group were concerned with HERA physics. In HERA 
electrons of E e w 27 GeV collide with protons of 
E p = 820 GeV , resulting in a centre of mass energy of 
y/s « 300 GeV . The kinematic region covered with the 
present data is roughly 10 -4 < x < 10 _1 , 7 GeV 2 < 
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Figure 1. Diagrams for DIS in 0(a a °) (quark partem model - 
QPM) and in 0(a s 1 ): boson-gluon fusion (BGF) and QCD 
Compton (QCDC) processes. 

Q 2 < 5000 GeV 2 and 40 GeV < W < 300 GeV , 
where W is the invariant mass of the hadronic system, 
Q 2 the negative 4-momentum transfer squared, and 
x the Bjorken scaling variable, to be identified with the 
proton momentum fraction carried by the scattering 
parton. HERA offers the opportunity to study the 
evolution of physics quantities over a large kinematic 
range. The large phase space available for hard QCD 
radiation, which can be treated in perturbative QCD, 
leads to prominent jets observable in the final state. 
Another area of recent interest is the kinematic regime 
at small x (x & 10~ 3 ) but sizeable Q 2 , not accessible 
at pre-HERA DIS experiments, where novel QCD 
dynamics are expected to play a role (e.g. pjj). 

The two HERA detectors ZEUS and HI |gj are large 
multipurpose, "almost 4-7r" detectors built around the 
beam line. Inner tracking detectors for charged particle 
detection are surrounded by a magnet and calorimetry. 
Both the scattered electron serving as a tag for DIS 
events and the hadronic final state are measured. Note 
that a substantial part of the hadronic final state, the 
proton remnant, leaves the detectors unobserved in 
the beam pipe. The region close to the proton beam 
direction is often referred to as the forward region. 

Apart from the laboratory frame, the hadronic 
centre of mass system (CMS) and the Breit frame are 
used in the analyses. The Breit frame is defined by 
the condition that the virtual photon does not transfer 
energy, only momentum. In the QPM picture the 
scattering quark would thus just reverse its momentum 
of magnitude Q/2. The CMS is defined as the centre 
of mass system of the incoming proton and the virtual 
boson, i.e. the CMS of the hadronic final state with 



invariant mass W. In both systems the hemisphere 
defined by the virtual photon direction is referred to 
as the current region, the other (containing the proton 
remnant) as the target region. The CMS current 
and target systems are back to back with momentum 
W/2 each. Longitudinal and transverse quantities 
are calculated w.r.t. the boson direction. With a 
longitudinal boost from the Breit frame into the CMS, 
particles formerly assigned to the target hemisphere may 
now end up in the current hemisphere. 

Monte Carlo (MC) models based upon QCD 
phenomenology are used to simulate the DIS process. 
The MEPS model (Matrix Element plus Parton 
Shower), an option of the LEPTO generator ||, 
incorporates the QCD matrix elements up to first 
order, with additional soft emissions generated by 
adding leading log parton showers. In the colour 
dipole model (CDM) [Q [f| radiation stems from a 
chain of independently radiating dipoles formed by 
the colour charges. Both programs use the Lund 
string model [|| for hadronizing the partonic final state. 
Deficiencies of the Herwig parton shower model have 
now been fixed by adding matrix element corrections 

[|. This model implements an alternative (cluster) 
fragmentation scheme ]l0| , allowing for valuable cross 
checks in the future. 

2. Jet physics 

The processes contributing to DIS up to first order in 
a s are shown in figure [| The QPM process results in a 
so-called "1+1" jet topology, while the QCDC and BGF 
processes give "2+1" jet events, where the "+1" refers 
to the unobserved remnant jet. From a measurement of 
the 2+1 jet rate at large x and Q 2 , where the parton 
densities are well known, a s can be measured. At small 
x and Q 2 , one can determine the largely unknown 
gluon density from the rate of 2+1 jet events, which 
is then dominated by the BGF graph (assuming a s to 
be known). Complications arise from the fact that 
the initial state contains strongly interacting particles, 
leading to the evolution of parton showers. Such effects 
need to be taken into account with the help of MC 
simulations. 

2.1. The strong coupling constant a s 

Both HI and ZEUS use the modified JADE algorithm 
[ pr[ with resolution parameter y cut — 0.02 to define 
jets in the a s analysis. A pseudoparticle is introduced 
to account for the unobserved remnant, and then all 
particles i,j satisfying mfj < y cut ■ W 2 are merged into 
jets. The chosen y cut value is a compromise between 
statistical precision (small y cu t ), and controllable higher 
order corrections (large y cut ) |12| . In the HI analysis 
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|l2) ] an angular cut 8j Ct > 10° (w.r.t. the proton 
direction) protects against parton showers close to the 
remnant. The obtained jet rates are corrected for 
detector effects, remaining parton shower contributions 
and hadronization with the MEPS model. In order to 
extract a s from the measured jet rates, it is important 
to take next to leading order (NLO) corrections into 
account to reduce dependencies upon y cut and the 
chosen renormalization and factorization scales 
Using PROJET @ as NLO calculation, the measured 
jet rate then yields measurements of a 3 (Q 2 ) in the 
range 10 GeV 2 < Q 2 < 3000 GeV 2 , which can be 
seen to run according to the QCD expectation fl3fl . 
However, below Q 2 = 100 GeV 2 , the corrections are 
very model dependent (MEPS vs. CDM). Therefore 
only data at Q 2 > 100 GeV 2 are used to extract 
a s (m z 2 ) = 0.123 ±0.018 g|. 

For 2+1 jet events with Q 2 > 160 GeV 2 and 
x > 0.01, ZEUS has measured the jet distribution in the 
Lorentz invariant z p variable ft5| , which in the centre 
of mass frame of the virtual photon and the incoming 
parton is an angular variable z p = \ ■ (1 — cos#j e t)- 
Here 0j e t is the angle of the jet w.r.t. the direction 
of the incoming parton. Perturbation theory in next 
to leading order (NLO) |l4|] is able to describe the jet 
angular distribution down to z p w 0.1. For z p < 0.1 an 
excess of jets is observed. Both, the MEPS (LO matrix 
element + parton showers) and ME (pure LO matrix 
element) simulations are similar to the NLO calculation 
|lq ]. The excess of jets at z p < 0.1 is therefore unlikely 
to be cured by next to NLO calculations. 

For the a s extraction, a cut z p > 0.1 restricts the 
data to a region well described by NLO perturbation 
theory and QCD models JTq| . The preliminary 
a s measurements @ for 100 GeV 2 < Q 2 < 
3600 GeV 2 demonstrate the potential of HERA to 
study the dependence of a s upon the renormalization 
scale, and agree well with the QCD expectation (see 
figure ||). It is expected that already the analysis of 
the 1994 HERA data, once finalized, will yield a very 
competitive measurement of a s (mz 2 ). 

2.2. The gluon density in the proton 

The 2+1 jet sample (defined with the cone algorithm in 
the CMS) in the range 10 GeV 2 < Q 2 < 100 GeV 2 
is used to extract the gluon density g(x g , Q 2 ) , because 
there the BGF graph (figure [l]) dominates (BGF:QCDC 
~ 4 : 1 jl6], [l?]])- The momentum fraction x g which the 
gluon carries is calculated from the invariant mass 2 s of 
the hard subsystem forming the 2 jets via x g — x(l + 
s/Q 2 ) « s/W 2 . Special cuts remove events affected by 
parton showers 0, 0. The MEPS model is used to 
unfold detector effects, the QCDC contribution, QPM 
background and remaining parton shower contributions. 
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Figure 2. Preliminary ct s (Q) measurements from ZEUS, 
compared to the QCD predictions corresponding to A^jg = 100, 
200 and 300 GeV. 

The MEPS model employs a cut-off for invariant parton- 
parton masses m% > y m i n • W 2 to regulate divergencies 
of its LO matrix element. In order to access x g as 
small as possible, s is chosen as small as experimental 
resolution allows, and as problems with the diverging 
LO matrix element can be avoided. It has to be 
ensured that the BGF events to be analyzed are actually 
generated by the model and do not fall below that cut- 
off @0. 

The HI analysis (l(| uses a fixed cut-off s > 
100 GeV 2 to define BGF events, and they parametrize 
the MEPS cut-off such as to follow the limit at 
which the order a s contribution exceeds the total cross 
section within a margin of AVI = 2 GeV . ZEUS 
uses the standard y m i n cut-off scheme in the MEPS 
model and defines BGF events via s > y m [ n ■ W 2 . 
The parameter y m ; n is then varied between 0.0025 and 
0.01 to study its influence on the result. The HI 
and ZEUS results |ll| agree well with each other, 
but yield different size systematic errors (figure ^). 
The ZEUS errors receive large contributions from 
the y m in variation. The rise of the measured gluon 
density towards small x can be described by a LO 
gluon density |l8| following the DGLAP (Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi) J2l| equations. The 
data are also consistent with the indirect determination 
of g(x g , Q 2 ) from the scaling violations of F 2 jlfj, 
providing a non-trivial test of QCD. 

2.3. Open Points 

Lack of understanding of parton showers close to the 
remnant (model dependent corrections, failure of NLO 
calculations) currently prevents the a s analysis to make 
full use of the large statistics data at Q 2 < 100 GeV 2 . 
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Figure 3. The gluon density in the proton, determined in 
leading order (LO) from the rate of 2+1 jet events. Shown are 
data from HI and ZEUS, compared to the LO GRV JJi) / gluon 
density parametrization. 
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Though increasing HERA luminosity will allow the 
a s analysis to be restricted to higher Q 2 to reduce 
uncertainties, the understanding of the forward region 
remains a challenge. 

So far the a s measurements rely solely upon the 
JADE algorithm, being the only algorithm for which 
NLO jet cross sections are calculated [|2]. NLO 
calculations for other algorithms, such as the cone [^3| 
or the theoretically preferred &t 0] algorithm are 
desirable. Such a program, which would also be able 
to calculate event shape variables like energy-energy 
correlations, Thrust, etc., is being worked upon by D. 
Graudenz, but results cannot be expected in a short 
term. Theoretical uncertainties could also be reduced 
by resumming higher order corrections. 

The validity of corrections from hadronic to partonic 
final states, defined either in LO or NLO, need to 
be checked with models based upon different parton 
shower and hadronization schemes. Unfortunately, a 
MC generator incorporating the QCD matrix elements 
beyond LO is missing. 

The gluon density has so far been determined in LO. 
A method allowing a measurement in NLO is presently 
under study Q. 

How can a s be determined consistently, considering 
it is input for the evolution of parton densities which are 
used in the analysis (25)? 

3. Novel QCD dynamics 

The observed strong rise of the structure function F2 
towards small x |^(| has caused much debate on whether 
the QCD evolution of the parton densities can still be 
described by the conventional DGLAP |2]J] equations, 
or whether the HERA data extend into a new regime 
at small x where the dynamics is governed by the 
BFKL (Balitsky-Fadin-Kuraev-Lipatov) |^6| equation. 
It would be extremely interesting to test QCD in such 



Figure 4. Parton evolution in the ladder approximation. The 
selection of forward jets in DIS events is illustrated. 

a new regime. While the rise is consistent with the 
expectation from BFKL dynamics, it can however also 
be described by a DGLAP evolution J27J. At lowest 
order the BFKL and DGLAP equations resum the 
leading logarithmic (a s hxl/x) n or (a s ln(Q 2 /Q 2 ,))™ 
contributions respectively. In this approximation the 
leading diagrams are of the ladder type (figure ||) • The 
leading log DGLAP ansatz corresponds to a strong 
ordering of the transverse momenta kx (w.r.t. the 
proton beam) in the parton cascade (Qq « ^ j « 
•••kx j <C ...Q 2 ), while there is no such ordering in 
the BFKL ansatz (kr | ~ hr [ p8| . Measurements 
on the hadronic final state emerging from the cascade 
therefore offer another handle to search for signatures of 
the BFKL behaviour. They are compared to analytical 
calculations as well as to the QCD models MEPS and 
CDM. The CDM description of gluon emission is similar 
to that of the BFKL evolution, because the gluons 
emitted by the dipoles do not obey strong ordering 
in k T The MEPS model with its leading log 

parton shower is based upon DGLAP dynamics, and 
the emitted partons are thus ordered in fey . 

3.1. Transverse Energy Production 

As a consequence of the strong kx ordering the DGLAP 
evolution is expected to produce less transverse energy 
Et in a region between the current region and the 
proton remnant than the BFKL evolution J3(|. HI and 
ZEUS have measured the flow of transverse energy in 
the laboratory frame as a function of pseudorapidity 
i] = — lntan(0/2), where 9 is the angle of the en erg y 
deposition w.r.t the proton beam axis [52J |33| . 
The measurements are made for varying ranges in 
x (2-10" 4 < (x > < 5-l(T 3 ) and Q 2 (7 GeV 2 < (Q 2 ) < 
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Figure 5. Transwr.se energy flow in the forward region at HI 
ply and ZEUS fell/ for x < 10~ 3 . The proton direction is to the 
right. The calorimeter acceptances end at n around 3.5. The 
data are compared to the CDM (here labelled CDMBGF) and 
MEPS models and to partonic calculations based upon the 
DGLAP and BFKL equations ml. 
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Figure 6. Transverse energy (Et ) per unit of pseudorapidity 
77* as a function of x for three different values of Q 2 , measured 
centrally at T)* = in the CMS ( corresponding to the lab. 
forward region). The data are compared to the CDM and MEPS 
models including hadronization, and to the BFKL calculation 
(no hadronization). 



30 GeV 2 and agree well between the experiments [ p2| . 

The Et flows for large x and Q 2 are reasonably 
well described by MEPS and CDM. For smaller 
x and Q 2 both models predict a more pronounced 
enhancement in the current fragmentation region than 
is seen in the data. Between the current system and 
the proton remnant (the lab. forward region), the data 
are reasonably well described by the CDM, while the 
MEPS model produces too little E T M. This 
intermediate region is expanded in figure |g, because 
there perturbative calculations, based either on DGLAP 
or on BFKL dynamics, are available f30[] . The BFKL 
calculation comes out close to the data, while the 
DGLAP calculation predicts much less Et ■ However, 
the non-perturbative hadronization phase is missing in 
these calculations. 

HI has determined the average Et , measured cen- 
trally in the CMS as a function of x and Q 2 (figure ^). 
They find an increase of (Et ) with decreasing x , which 
is a characteristic BFKL prediction p0| . The data are 
in agreement with the BFKL calculation (3^], if one as- 
sumes an Et contribution from hadronization of about 
0.4 GeV per unit rapidity (independent of x ). That 
estimate is taken from the CDM, which agrees with the 
BFKL calculation at the parton level. 

The apparent failure of the MEPS model has caused 
many questions about its ingredients: the way the 
parton shower is "matched" to the matrix element, 
the colour connection between the current and the 
remnant system and its effect upon hadronization, 
and the remnant fragmentation itself which is little 
tested. It seems that re-arranging colour configurations 



can produce enough Et through hadronization to 
compensate the Et deficit in the DGLAP cascade 
of the MEPS model |§. A MEPS version thus 
modified should be available soon for detailed testing. 
The flexibility in the hadronization modelling presently 
precludes unambiguous tests of the DGLAP evolution 
through Et measurements. For the same reasons the 
intriguing success of the CDM without Ut ordering may 
be fortuitous. A MC model invoking explicitly the 
BFKL evolution, currently being developed by K. Golec- 
Biernat et al., would help interpreting the data. In 
any case, the Et data provide important input for QCD 
phenomenology. 

3.2. Forward Jets 

At present strong conclusions upon the validity of the 
BFKL or DGLAP parton evolutions at small x from the 
Et measurements are hampered by the uncertainties 
about hadronization. Jet production should be less 
affected by hadronization. A signature for BFKL 
dynamics proposed by p6[ is the production of "forward 
jets" with Xj c t = Ejct/Ep, the ratio of jet energy and 
proton beam energy, as large as possible, and with 
transverse momentum fcrjet close to Q in order to reduce 
the phase space for the /ct ordered DGLAP evolution 
(see figure ||). An enhanced rate of events with such 
jets is thus expected in the BFKL scheme |}6|, [3?J . The 
experimental difficulty is to detect these "forward" jets 
which are close to the beam hole in proton direction. 

The rate of forward jets measured by HI ||| |3^] 
(figure 0) is larger at low x than at high x. This 
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Figure 8. The laboratory angular distribution of jets detected 
either in the Breit current or target hemisphere. 



Figure 7. The rate of forward jets (selected with xj ct > 0.025, 
0.5 < fcTjet 2 /Q 2 < 4 and fcyjet > 5 GeV ) in the kinematic 
range 2-10" 4 < x < 2-10" 3 and Q 2 » 20 GeV 2 . The 
measurement is compared to the CDM and MEPS models. 



jets. If these jets can be identified with gluons emitted 
from the ladder, it would be possible to check the parton 
ordering directly. 



is expected from BFKL calculations, in contrast to 
calculations without BFKL ladder |37|, The 
behaviour of the data is better represented by the 
CDM than by the MEPS model. However, neither of 
them describe the energy spectrum of the observed jets 
correctly, and the model pre dictions for the jet rates are 
thus cut dependent ]38|, p3| . The analysis of a larger 
statistics sample should allow more firm conclusions. 

ZEUS has measured an inclusive jet cross section 
dcr/d-Erjet in the Breit frame |3q| . Many more jets are 
found in the target region with a harder -Erjet spectrum 
than in the current region, reflecting the differences in 
phase space in the two regions. The current region 
data are reasonably well described by the CDM and 
MEPS models. In the target region however there is 
a substantial excess of jets over the model predictions, 
which can be linked with an excess of 2 + 1 jet events 
|3§| ]. In the laboratory frame this excess is located in 
the forward region at angles 9 jct < 20° (figure ||) . 

The data on jet production in the forward region 
(lab. frame), or the target region (Breit frame) certainly 
pose a challenge to theory. So far cross sections are 
calculated J3l], [39| only for partons, while experiments 
measure hadron jets. This gap has to be bridged from 
both sides to allow a strictly valid comparison. 

3.3. Jet correlations 

Apart from calculations of forward jet rates, Del 
Duca [[39| discussed angular correlations for forward 
jets. If a BFKL ladder is inserted in between the 
electron-photon vertex and the forward jet, the angular 
correlation between the forward jet and the electron 
imposed by momentum conservation is relaxed. Such 
angular decorrelation could be another footprint of 
BFKL dynamics. The fact that 4% of the HI forward jet 
events contain a second forward jet |3S|l opens another 
route of investigation, namely correlations between such 



3.4- Dipole emission 

An interesting ansatz to calculate final state observables 
was presented by R. Peschanski p2|. The starting 



point is onium-onium scattering 43 with onium wave 
functions which can be derived from QCD. Such a 
reaction is analogous to an interaction of the current 
system with the remnant system in DIS. Radiation 
is treated in the dipole picture, leading to a copious 
production of dipoles in the central rapidity region of 
the interaction. Once such an ansatz yields quantitative 
predictions, it could be tested in DIS, e.g. with E T flow 
measurements. 

Bo Andersson jl0| discussed DIS final states in terms 
of a chain of radiating colour dipoles, and its connection 
with the Ciafaloni-Catani-Fiorini-Marchesini ansatz 
[|4~i| . In principle this model could provide a complete 
picture of the hadronic final state in DIS. The 
implementation in the Ariadne || MC generator is in 
progress to allow detailed predictions. 

3.5. QCD Instantons 

The standard model contains processes which cannot 
be described by perturbation theory, and which violate 
classical conservation laws like baryon and lepton 
number in the case of the electroweak sector and 
chirality for the strong interaction [ }44j . Such anomalous 
processes are induced by instantons jf5|. At HERA, 
QCD instantons may lead to observable effects in 
the hadronic final state in DIS |47]], which were 
discussed by F. Schrempp. The instanton should 
decay isotropically into a high multiplicity state of 
gluons and all quark flavours simultaneously which are 
kinematically allowed. A MC program to simulate 
instanton events has become available Due to the 
isotropic decay, one expects a densely populated region 
in rapidity, other than the current jet, which is isotropic 
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in azimuth. The presence of strangeness and charm 
could provide an additional signature. However remote 
the a priori chances to see such signals may appear, here 
is a chance for a major discovery at HERA! 

4. Charged Particle Spectra 

The HI and ZEUS measurements of inclusive charged 
particle spectra Q are performed either in the Breit 
frame or in the CMS. In the Breit frame in- and outgoing 
quark have equal but opposite sign momenta Q/2 (QPM 
picture), and in e + e~ annihilation the outgoing quark 
and antiquark have equal but opposite momenta yfs/2 = 
Q/2. Due to this similarity it is interesting to compare 
particle spectra in the Breit current hemisphere in DIS 
with e + e _ data. DIS experiments have the advantage 
over e + e~experiments that they cover a large span in 
Q , presently from 3 GeV to 50 GeV , in a single 
experiment. The current mean charged multiplicity 
at HERA rises ~ InQ within errors, and agrees with 
e + e~data (divided by 2) where they overlap [[50[ [49| . 

Colour coherence should lead to a suppression of soft 
gluon emission. The HERA data j50[ on the scaled 
charged particle momentum distribution In 1 / x p with 
x p = 2 ■ p/Q exhibit the expected hump backed plateau 
|M , the evolution of which with Q is in agreement with 
the assumption of colour coherence. However, like in 
e + e~annihilation, this behaviour can also be mimicked 
through the Lund string fragmentation ]49[ . 

The scaled momentum spectrum of Xp in the CMS, 
where the particle longitudinal momenta p z are divided 
by the maximal possible momentum, xp — 2 -p z /W, 
are shown in figure ^| for the current region (the target 
region is not observed). Comparing HERA data at 
W w 120 GeV M w ith fixed target data at W = 
14 and 18 GeV |53, [52f, significant scaling violations 
are observed, in agreement with QCD expectations: the 
large value of W at HERA results in a large phase space 
for QCD radiation, softening the xp spectrum w.r.t. 
data at lower W . It can be expected that such data 
will be used to extract a s in the future. 

The effect of QCD radiation is clearly seen in the 
"seagull plot" (figure |lC| ), where the mean transverse 
momenta p\ squared of the particles is plotted as a 
function of xp . As a consequence of increased QCD 
radiation, much larger p\ are observed at HERA |jl], [l9| 
than at EMC |g3| at smaller W, again in agreement 
with QCD expectation. ZEUS has also compared DIS 
events with and without a large rapidity gap [ pi[ in 
this respect 49 . Much smaller p\ than in normal 



DIS events are observed in events with a large rapidity 
gap, thought to stem from diffractive processes and 
accounting for approximately 10% of the total sample 
|H^1 . This indicates that the scale governing radiation is 
much smaller than W for rapidity gap events. 
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Figure 9. The xp spectra measured at HERA compared with 
the QPM (dotted line) without QCD radiation, the MEPS model 
(full line), and with fixed target DIS data at lower W . 



5. Conclusion 

Two complementary approaches to the HERA data can 
be distinguished. In one approach, one tries to identify 
a region which is "well understood" , meaning that the 
observation agrees with the theory and the models. 
Under this condition, the data can be interpreted in 
the framework of the theory, and physical quantities 
which are defined within the theory can be extracted. 
The measurements of a s and g{x g ,Q 2 ) fall into this 
category. However, we have also seen data which are 
not yet understood theoretically, namely hadron and jet 
production in the forward region. Such data currently 
pose a challenge to the theory, and experimentalists 
should make every effort to provide theory with solid 
data to work with. 
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HERA data with and without a rapidity gap (LRG/NRG) 
compared to the QPM prediction (dotted line) and the MEPS 
model (full line), and to EMC data at lower W . 
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